The cool end of the DZ sequence in the SDSS 
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Abstract. We report the discovery of cool DZ white dwarfs, which lie in the SDSS (u-g) vs. (g-r) 
two-color diagram across and below the main sequence. These stars represent the extension of the 
well-known DZ sequence towards cooler temperatures. 
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INTRODUCTION 

Cool He-rich metal polluted white dwarfs (spectral class DZ) in the SDSS data have 
been analyzed by [1]. In the (u-g) vs. (g-r) diagram (Fig. 1) all but one fall in the region 
above the main sequence of normal stars (grey dots). Of the 147 objects discovered in 
the SDSS only two have T eff below 6600 K (at 6090 and 4660 K). This suggests that 
a large number of cooler DZ remains to be found. We have therefore visually checked 
all spectra in the region of QSO colors with borders described by the dashed lines and 
identified 17 as cool DZ. The location in the two-color diagram indicates that the DZ 
sequence might "tunnel" through the main sequence. Using a new search routine in the 
region bounded by the continuous lines, which checks for the characteristic features 
(most prominently Mgl lines near 5 170 A), we found another 6 objects, 2 of which were 
already in the Dufour sample. 

Fig. 2 shows a selection of typical spectra. SDSS 1535+1247 (= NLTT40607 = G137- 
24) was already identified as a DZ and analyzed by [2] and [3]. It looks extremely 
similar to G165-7, shown by [4] to be weakly magnetic. That star is also in our sam- 
ple (SDSS1330+3029). SDSS0916+2540 has the strongest features of any DZ known, 
whereas the remaining two show little more than the strong asymmetric feature between 
5000 and 5200 A, which is characteristic for this group. SDSS 1336+3547 is magnetic, 
with multiple line cores in many metal lines. 

Our list of extreme DZ is given here 

SDSS 0157+0033 SDSS 0205+2155 SDSS 0916+2540 SDSS 0925+3130 

SDSS 0937+5228 SDSS 0956+5912 SDSS 1033+1809 SDSS 1038-0036 

SDSS 1040+2407 SDSS 1043+3516 SDSS 1103+4144 SDSS 1144+1218 

SDSS 1152+1605 SDSS 1234+5208 SDSS 1330+3029 SDSS 1336+3547 

SDSS 1404+3620 SDSS 1421+1843 SDSS 1430-0151 SDSS 1524+4049 

SDSS 1535+1247 SDSS 1546+3009 SDSS 1616+3303 
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FIGURE 1. SDSS Two-color diagram for objects with spectra identified as normal stars (grey S-shaped 
region in the middle of the figure), QSO (dots below main sequence), DZ white dwarfs analyzed by [1] 
(heavy black dots) and cool "extreme" DZ discovered by our search (small circles). 

ANALYSIS 

All spectral features can be identified with lines from Ca, Mg, Na, Fe, Ti, and Cr. 
The broadest lines show strongly asymmetric profiles, which were in the case of Mgl 
5169/5174/5185 correctly identified as due to quasistatic broadening by [2]. The width 
of these lines, as well as that of the Cal and Call resonance lines is far beyond the range 
of validity of the impact approximation. We have used the simple and elegant method of 
[5], who present numerical calculations for the transition range between impact and 
quasistatic regime, which can in both limits be easily extended with the asymptotic 
formulae. These profiles are reasonable approximations in many cases. Nevertheless, 
the current analysis is only a first attempt and the results are a preliminary estimate of 
the general parameters of these objects. Among the problems which will need to be 
resolved are 

• The Cal and Call resonance lines between 3900 and 4250 A are so strong that 
they completely cut off all the flux below 4500 A. Our profiles assume a van der 
Waals r~ 6 dependence for the perturbation energy, which has been shown by [6] 
to be a very poor approximation for the Ca-He interaction. This assumption has to 
be replaced by a more sophisticated potential for each individual line. The current 
profiles are not a satisfactory approximation. 

• The u-g colors suggest that the UV flux near 3500 A is almost totally suppressed in 
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FIGURE 2. Four stars from our sample of "extreme" DZ stars (blue), together with theoretical models 
(red), which at least qualitatively reproduce the line features and the continuum slope. If the models are 
not visible in black and white print, they fall exactly on the observed spectrum. 



all objects. We have included all expected strong lines in the UV - in particular from 
Fel and Mgl/Mgll - in our modeling, which results in a strong blanketing effect 
in the optical. The models qualitatively reproduce the u-g colors, and therefore 
also the blanketing effect. However, the Mgl/II resoncance lines are expected to 
be extremely strong, and similar caution applies to their profiles as for the Cal/II 
resonance lines. 

The upper limit to the H/He ratio is ps 10 4 from the non-visibility of Ha, but the 
exact amount is unknown. Since the broadening of Lyman a by helium leads to 
very strong absorption even into the optical range, this is a further unknown free 
parameter. 



• Finally, the surface gravity can also not be determined from the available observa- 
tions, and we have assumed logg = 8. 

RESULTS 

For a preliminary analysis we calculated a grid of helium-rich atmosphere models with 
T e ff between 7000 and 5000 K, log g fixed at 8.0, and variable metal abundances 
of 13 most commonly detected elements between Na and Fe, with the ratios kept at 
approximately solar values. The hydrogen abundance was fixed at log H/He = -4. A 
temperature was determined from the slope of the spectrum, which is in most cases 
confirmed by the SDSS magnitudes. In some cases additionally the ionization balance 
of Cal/Call was taken into account. The objects in the two-color diagram below the 
main sequence are essentially the continuation of the DZ temperature sequence down to 
« 5400 K. 

For the most metal-rich objects individual abundances were adjusted until a reason- 
able fit to the spectral features was achieved. As an example we note results for two 
objects with the strongest features 

SDSS 1535+1247: T eff = 5700 K, [Mg/He] = -7.50, [Na/He] = -8.90, [Ca/He] = -8.90, 
[Fe/He] = -7.60 

SDSS 0916+2540: T eff = 5500 K, [Mg/He] = -6.90, [Na/He] = -9.25, [Ca/He] = -7.50, 
[Fe/He] = -7.10. 

It is obvious that - while the Mg/Fe ratios are quite similar - the relative Ca and Na 
abundances vary significantly. Since the diffusion time scales of these elements agree 
within about 30%, it is quite likely that these differences have to be attributed to different 
compositions of the accreted material. While the Ca abundance of SDSS 0916+2540 
looks exceptionally high compared to the Dufour et al. results, the other objects follow 
the general trend of their Fig. 9. 

The total mass of observed metals in the convection zone of SDSS 0916+2540 is 
~ 5.2 10 21 g. Adding C, O, Si with the relative abundances as observed in GD40 [7] the 
total mass of metals is ~ 1 .4 10 23 g. If accretion really has reached diffusion equilibrium 
(the diffusion time scales are typically 300000 yrs) the total accreted mass is of the order 
of the Pallas or Vesta masses. 
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